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ABSTRACT 



A process and related apparatus for producing at least two 
matched color displays of a digital image using two different 
display devices, comprising using an adaptor to convert the 
digital information representing the image to digital infor- 
mation such that the displayed image as a result of this 
converted digital information on one of the devices, appears 
the same as the image displayed on the other to an average 
observer when viewed under similar viewing conditions. 
The adaptor contains a transform LUT constructed using 
display device modeling to generate corresponding values to 
preselected input digital image values to map the input color 
values into the color space of the display device so that 
image matching is achieved. A output value signals received 
from a first model are compared to output value signals 
received from a second model to obtain an error signal used 
to modify the input color values to the second model, and the 
modified color values are then used to create die transform 
LUT for matching the color display of the second display 
device with that of the first device. 



12 Claims, 6 Drawing Sheets 
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COLOR MATCHING METHOD AND 
APPARATUS 

This is a continuation of application Ser. No. 07/992,874 
filed Dec. 15, 1992, now abandoned. 5 

COPYRIGHT NOTICE 

A ponion Of the disclosure of this patent document 
contains material which is subject to copyright protection. 10 
The copyright owner has no objection to the facsimile 
reproduction by anyone of the patent document or the patent 
disclosure, as it appears in the Patent and Trademark Office 
patent file or records, but otherwise reserves all copyright 
rights whatsoever, 15 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

This invention relates to color matching and more par- ^0 
ticularly to a process and related apparatus for matching the 
color displayed on a plurality of color display devices. 

2. Description of Related Art 

The primary use of this invention is in the field of 
computer assisted color publishing systems particularly in 
the area of color matching. In such systems typically a color 
image is scanned using a scanning device which measures 
light intensity reflected or transmitted from the image to 
generate electronic signals representing three achromatic 
images of the original, each of the three achromatic images 
resulting from detecting the light intensity from the original 
color image through a colored filter. The filters used are 
almost always Red, Green and Blue. The electronic signal 
from the scanner is converted to a digital signal in which 
light intensity levels are represented as numbers. Informa- 
tion identifying each set of numbers representing the image 
information obtained through each filter is also preserved. 
Thus, through these steps, the original colored image is 
converted to a plurality of image values, and for each picture 
elemenl in the picture, there are three such values, a Red, a 
Green and a Blue. 

The image represented by the image values may be 
displayed in a CRT type monitor, or may be printed using a 
printing device able to accept electronic input. Often the 45 
image will be displayed at different locations and times 
using more than one display device, such as a plurality of 
CRT displays for observation and study by more than one 
operator. Similarly, hard copies may be desired in two or 
more diflFerent locations using two or more different printers. 50 
Even though the input values to the multiple CRT displays 
are the same, the displayed colors are vastly different, as 
anyone is well aware of who has ever observed the multiple 
television displays in a typical television sales store front. 
However, when one needs to make decisions about accept- 55 
ability of color for a display, one needs to know with a great 
degree of confidence that the color, that one is observing and 
discussing with an operator who is observing the same 
image on a different CRT, is the same as the color observed 
by the operator. The same is true, if the displayed image is 50 
one created on a printer and compared with the same image 
created on a different printer. 

Colorimetry, which is the study of color based on both 
spectral disuibution of the energy reflected or transmitted 
from a test sample and the response of the human eye, as 65 
well as the special distribution of the illuminating source, 
provides a method to describe and measure color and 



enables one to determine when colors match. Through the 
use of CIE defined TCstimulus Values (which are the 
amounts of three primary lights which when added produce 
a visual, or colorimetric match with an original color), one 
may determine with reasonable certainty that if two colors 
have the same three CIETristimulus values that is, if the Red 
Tristimulus value of one color is the same as the Red 
Tristimulus value of the other, and so on for the Green and 
Blue Tristimulus values, then the appearance of the two 
colors will be the same to the average observer. CIE stands 
for the International Commission on Dlumination. 

It appears, therefore, that the problem of matching color 
outputs of different displays is readily resolved by providing 
displays having the same CIE tristimulus values for the same 
image value inputs. 

While the solution in principle is simple, creating a 
conversion device which will make two or more different 
displays produce the same tristimulus output values for the 
same input image values is very difficult. Each of the display 
devices operates in its own, device dependent, color space 
where image values at its input are transformed into display 
image values. The difiSculty lies primarily in matching the 
two transformations occurring within the two display 
devices for the image values appearing at the input of each, 
so that both display the same displayed colorimetric values 
for the same input image values. 

The prior art solutions to color matching fall into two 
fundamentally distinct approaches. The first is based on the 
decomposition of a color vector to a set of primaries, and 
known as the primary decomposition technique. Typical of 
this approach is the technique described in U.S. patent 
application 5,196,927 filed Mar. 16, 1990, assigned to E. I. 
du Pont de Nemours and Company. According to this 
technique, the input image is decomposed into the linear 
combination of eight primaries (yellow, magenta, cyan, 
black, red. green, blue, and three-color black). CoefBcients 
for a 4x8matrix (the values of CYMK for each of the eight 
primaries) are adjusted such that the color of the eight 
primaries in the input from both systems to be matched are 
matched, and the new input values are found by matrix 
multiplication of an 8x1 matrix (coefficient of decomposi- 
tion) with the 4x8 matrix. This technique suffers because of 
the non linearity in the additive properties of color dyes and 
because the "primaries" used are not mathematically inde- 
pendent. 

The more common solution, adopted by the printing 
industry is the grid sampling technique. This involves using 
a color transformation formula such that the error between 
the targeted color and the processed color is minimized. The 
differences between the various methods in existence are 
found in the specific transformation formulas. U.S. Pat, No. 
4,500,919 issued to Schreiber is a good example of the use 
of a transformation formula. The difficulty with this 
approach again lies in the non-linearity of the color addition 
process and the complexity of the color surfaces in a set of 
equations. 

SUMMARY OF THE INVENTION 

The present invaition relates to provide a method for 
matching the color display of at least a first and a second 
display devices, comprising: 
(I) creating a transform LUT for converting input color 
values to output color values by: 
(1) producing a first preselected plurality of input color 
values; 
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(2) using said plurality of input values to display a 
plurality of color patches in said ftrst and said second 
display devices; 

(3) obtaining a colorimetric value of each of the dis- 
played color patches in each of the displayed devices 5 
and using said colorimetric values to create a first 
and a second model for the first and second devices 
respectively correlating preselected input color val- 
ues to displayed colorimetric values for each of the 
two devices; jq 

(4) inputting to said first and second models a second 
preselected plurality of color values; 

(5) comparing the output of the first model to the output 
of the second model to obtain an error signal indica- 
tive of the difference between the two output signals; 

(6) using the error signal to modify the input color ^5 
values to the second model and again comparing the 
output of the first model to the output of the second 
model to obtain a new error signal; 

(7) repeating the process of steps (5) and (6) above until 
the error signal is a minimum; and 20 

(8) using the modified color values to create a trans- 
form LUT correlating input values to modified val- 
ues; 

(II) using the transform LUT to transform any plurality of 
source color image values before inputting said image 25 
values to an input of the second display device; and 

(ni) displaying said transformed source color image val- 
ues on said second display device. 

When there is no corresponding source image value in the 
transform LUT correlating an input source color image 30 
value to a modified color image value, one can use inter- 
polation to derive a corrected source image value from a 
closest source and corresponding modified value in the LUT. 

Preferably, the first and second preselected pluralities of 
color image values are the same. 35 

This invention further relates to an apparatus for gener- 
ating a transform LUT for converting input color values to 
output color values comprising: 

(1) digital color image values input means; 

(2) first and second display model means for producing 40 
each an output colorimetric tristimulus image values 
for input color values, each of said model means 
having an input and an output, the second model 
input connected to the input means; 

(3) an adder device having a first signal input coimected 45 
to the input means, a second, correction error signal 
input, a control signal input, and an output connected 

to the first model input, for outputting modified color 
values; 

(4) means connected to the output of said first and 50 
second model means, for comparing the tristimulus 
values output of the first and second display models 
and for producing an error signal; 

(5) means for testing the error signal to determine if 
said signal is a minimum error signal, and for out- 55 
putting a correction error signal and a control signal 

to the adder, said testing means cormected between 
said means to compare and said adder device; and 

(6) means also connected to the input means, the adder 
output, and the means for testing, for receiving the 60 
adder output and for generating a transform LUT 
correlating digital input color values to the modified 
color values appearing at the output of the adder 
when the means for testing determines that the error 
signal is a minimum. 65 

The apparatus may further comprise means for storing the 
generated LUT. A gamut mapping means may be included 



between the output of the second mode! and the means to 
compare for mapping output values from the second model 
into a color gamut commensurate to that of the first model. 

The apparatus may comprise hardware, or may be a 
computer programmed through software to perform all of 
the above operations, or may be a combination of dedicated 
hardware and computer implemented software. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention can be more fiilly understood from the 
following description thereof in connection with the accom- 
panying FIGURES described as follows. 

FIG. 1 is a schematic representation of an arrangement in 
which two matched images are displayed on two CRT 
displays. 

FIG. 2 is a schematic representation showing the arrange- 
ment used in obtaining data for generating the LUTs used in 
modeling the two displays of FIG. 1. 

FIG. 3 is a schematic representation of the process for 
obtaining the transform LUT used in the adaptor in the 
arrangement shown in FIG. 1. 

FIG. 4 is a generic representation of a display model. 

FIG. 5 is a schematic representation of the process used 
in deriving the correction factors shown in FIG. 3. 

FIG. 6 is a schematic representation of an arrangement in 
accordance with this invention in which two matched 
images are printed using two printing devices. 

FIG. 7 is a schematic representation showing the arrange- 
ment used in obtaining data for generating the LUTs used in 
modeling the two printing devices of FIG. 6. 

FIG. 8 is a schematic representation of the process for 
obtaining the transform LUT used in the adaptor in the 
arrangement shown in FIG. 6. 

• FIG. 9 is a schematic representation of an arrangement in 
accordance with this invention in which two matched 
images are printed using a CRT display and a printer 
including an RGB to CMYK converter. 

FIG. 10 is a schematic representation of an airangemem 
in accordance with this invention in which two matched 
images are printed using a CRT display and a color proofing 
device using a YMCK to RGB to drive the CRT display. 

FIG. 11 is a generic representation of a display model for 
use in developing the transfer LUTs required for the embodi- 
ments shown in FIGS. 9 and 10. 

DESCRIPTION OF THE PREFERRED 
EMBODIMENT(S) 

The invention will next be described in detail with refer- 
ence to the drawings in which similar characters indicate 
similar elements in all figures of the drawings. During the 
following discussion, color values are quantized in 256 steps 
corresponding to an 8 bit system. Other quantizations are 
possible, and not intended to be excluded because of the use 
of the 256 steps in the following examples. 

The invention comprises a method and apparatus for 
providing color visual matching of two representations of 
the same image when the image is displayed in two different 
displays, so that an observer will on visual examination of 
the displayed images observe two images having substan- 
tially the same color appearance. Using the method or 
apparatus of this invention causes each of two observers, 
whose visual response substantially conforms to the visual 
response of an average observer as defined by the CIE 
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institute, upon observing any one of the color displays under 
the same or similar surrounding illumination and back- 
ground, to receive substantially the same visual impression. 

FIG. 1 shows in schematic representation a situation 
where a color image is displayed on two CRT type displays 5 
which receive RGB type inputs to display an image. The 
image may be stored and manipulated in digital format in a 
work station 10 which itself may comprise a scanner, an 
image processor, a display device, and have image storage 
capabilities, or any portion of the above, or more. jq 

Because the required input to displays 12 and 14 in this 
example is a digital RGB type signal, the work station 
outputs over line 16 a digital RGB signal representing a 
color image. The digital RGB signal comprises a set of color 
values, one set of three values (R, G, B) for each picture 15 
element of the image 10 be displayed. The signal is directed 
over line 18 to the first display monitor 12 on which the 
colored image is displayed. The digital RGB signal is also 
sent over line 20 to the second display monitor 14, which 
could be located at a location different than the location of 20 
the first display 12. 

Prior to applying the RGB signal to the second monitor 
14, there is interposed an adaptor 22 which receives the 
digital RGB signal and after processing the signal, outputs 
a signal comprising a new set of digital color values R'. G', 25 
B' over line 24 to an input of the second display monitor 14 
for each image pixel. 

The adaptor 22 comprises a receiving and determining 
element 26 which receives each set of R, G, B values and 
checks them against a plurality of sets of R, G, B values 30 
comprising a transform look up table (LUT) 28 which 
correlates R, G, B values to a new set of R', G', B' values 
which when used as an input to the second display monitor 
14 will generate a pixel whose colorimetric value will match 
that of the pixel generated by the same R, G, B set of color 35 
values in the first display monitor 12. 

If the receiving and determining element 26 identifies a 
corresponding set of R, G, B values in the LUT 28, it outputs 
the appropriate R*, G', B' values on line 24. If not, the R, G, 
B input is directed to an interpolator 30 where R', G', B' ^ 
values arc derived by interpolation using existing adjacent 
sets of R, G, B values in the LUT, Preferably, the selection 
of the R» G, B values for the development of the transform 
LUT 28 is such that values derived by interpolation produce 
a visually acceptable match of the displayed pixel in the two 
monitors 12 and 14. 

For this system to operate with any degree of success, the 
development of the transform LUT 28 is critical. The values 
generated by this LUT 28 must indeed lead to values that 
correctly reproduce colors on the second display 14 which 
visually match the colors on the first display 12. FIGS. 2, 3 
and 4 help explain how the LUT 28 is generated. 

The first step in developing the transform LUT 28 is the 
creation of two models representing the two displays whose 55 
output is to be matched. FIG. 2 shows how this is accom- 
plished. 

The work station 10 is used to generate and output over 
line 16 a set of digital R. G, B color values in a regular 
pattern which is selected to supply color values closely 60 
enough spaced so that interpolation of values between 
existing values is reasonably accurate. Preferably, but not 
necessarily, the output is adjusted so that the color values 
produce displays 32 and 34 on the first and second display 
monitors 12 and 14, respectively, comprising a plurality of 65 
different color patches on display screens of the display 
monitors 12 and 14. As an example, the work station 10 may 



6 

generate a plurality of digital RGB values such as shown in 
the following table I: 

TABLE I 

R: 0, 13, 26, 51, 76, 102, 128, 153, 178, 204. 230, 255. 
G: 0, 13, 26, 51, 76. 102. 128, 153, 178. 204. 230, 255. 
B: 0, 13. 26, 51. 76, 102, 128. 153, 178, 204, 230, 255. 



Each combination of R, G, B values from this table 
represents a set (R, G, B)„ of color values. In this illusUation 
where there are 1728 possible such sets of R, G, B values, 
n=l to 1728. 

These 1728 sets of R. G. B values are supplied to both 
display monitors 12 and 14 over lines 18 and 20' creating a 
total of 1728 displayed patches for each of the monitors 12 
and 14, corresponding to the 1728 different combinations of 
the selected values. All 1728 patches are not necessarily 
displayed on each monitor screen simultaneously. 

A colorimetric measuring device 36 is used to read the 
displayed patches displayed in the first and second monitors 
12 and 14 as shown in FIG. 2. The output of the colorimetric 
device is a colorimetric set of color values for each of the 
patches. In our preferred embodiment, the output of the 
colorimeter are CIE defined L, a and b sets of color values 
in the Lab color space. Output colorimetric measurements 
given in other colorimetric color spaces are equally accept- 
able, such as XYZ, Tristimulus CIE defined RGB, etc. 
Uniform as well as non-uniform color spaces may be used. 
However, values in non-uniform color spaces related to a 
uniform color space with a known mathematical relation 
may result in increased calculation steps which must be 
compared with any gained advantage as a result of this 
selection, to determine if the particular choice of color space 
is justified, (For a description of the different color spaces 
and related terminology, see, in general, DIGITAL IMAGE 
PROCESSING, 2nd edition, by Mlliam Pratt, published by 
John Wey and Sons, Inc. pages 62-73.) 

Once the measurement of all patches is completed, two 
LUTs are compiled, representing the transfer functions of 
the two monitors 12 and 14. The first LUT will consist of the 
(RGB)„ sets of values and the corresponding (Lab)„ sets of 
values read off the display on monitor 12 and the second 
LUT will consist of the same (RGB)„ sets of values and the 
corresponding CL'a'b')„ sets of values read off the display 
monitor 14. 

Referring to FIG. 3, a first and a second model 40 and 42 
are built representing the two monitors 12 and 14, respec- 
tively. FIG. 4 shows a generic model 44 structure, which is 
used throughout this description in this invention to convert 
color values. The model 44 comprises an LUT 48, which is 
the particular device derived LUT representing the device 
transfer function. In this example, the LUT 48 for the first 
model 40 will be the first LUT correlating the (RGB)n sets 
of color values to the (Lab)„ sets of values, and the LUT 48 
for the second model 42 will be the second LUT correlating 
the same (RGB)„ sets of color values to the (L'a'b')n sets of 
values. 

In addition to the LUT 48, the model 44 comprises a 
receiving and determining means 46 similar to, or the same 
as, the receiving and determining means 26 described ear- 
lier, and a mathematical interpolating means 50, also similar 
to, or the same as, the previously described interpolating 
means 30. 

HG, 3 schematically represents the generation of the 
transform LUT. Sets of R. G. B color values from a 
preselected plurality of R, G, B sets of color values are used. 
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This plurality of R, G, B sets of values is, preferably, also 
produced in the work station 10, and, again preferably, is the 
same as the sets of (R, G, B)„ values from Table I previously 
used to generate the two LUTs for the two display models 40 
and 42. 5 

Each (R. G, B) n set of values is directed to the input of 
the first model 40 over lines 54 and 58. Model 40 produces 
an output of (Lab)n values corresponding to the input (R, G, 
B)„ values for this model 40. The same (R, G, B)„ values are 
directed to the second monitor model 42 over line 56. Ahead 
of model 42, there is an adder 60 which operates to add to 
the R„, G„, and B„ components of the (R, G, B)„ set of values 
any correction factor dR„, dG„, dB„, appearing on line 78. At 
first, nothing appears on line 78, and the R„, G„ and B„ 
values are appUed to the input of model 42 unaltered. Model 15 
42 also produces an output set of (L'a'b')„ color values 
corresponding to the input (R, G, B)„. This output appears on 
line 64. 

The (Lab),, values over line 72 and the (L'a'b')„ values 
over line 64 are next compared in comparator 66 and a 
difference signal [(Lab)„-(L'a'b') J is generated and directed 
over line 74 to correction factor generator 76, Using this 
difference signal, correction factor generator 76 produces 
correction factors dR^^j^, dG„(ij, ^^no) (^'^ ^ manner to be 
explained herein below) which arc sent over line 78 and ^ 
added to the corresponding components of the (R, G, B)„ set 
of values in adder 60, resulting in a new set of color values 
Rnci;=R«+dRnci>. G„^i;=G„+dG,(ij, and B„(,;=B„-KlB„(jj 
on line 62. These new values are again applied to the input 
of model 42 which produces a new output set of Lj'aj'b^ 
values on line 64 which is again compared in comparator 66 
with the set of (Lab)„ values from model 40. A new set of 
correction factors dR„(2). ^^ni2y ^^d dB„(2) is produced and 
added to R„(i/, G„(i)' and B„(,j' to produce a new set of R„', 
G„', and values An error "E" defined as: 



30 



35 



E=[ a L')V(a— a')^+{l>— b')21 

is used to determine when to stop this cycle. Referring to 
FIG, 5, error "E" is compared in comparator 92 with 40 
preselected minimum acceptable limits or a "0'* value. If the 
error is "0", or within the preselected acceptable limits, a 
signal over line 94 switches switch 98 (which is part of adder 
circuit 60) to feed the R„', G„,' and B„' values which 
produces the minimum over line 80 to the transform LUT 45 
compiler 82, rather than to model 42; the next input set of 
values (R, G, B)^^ is then applied over line 54 to the two 
models and the whole process is repeated for this new set of 
values and so on. 

The set of R„'G„'B„' values which produced this zero or 50 
acceptable minimum error are sent over line 80 to the LUT 
compiler 82 where the (R, G, B)„ and corresponding 
(R'G'B')„ sets of values are used to compile the transform 
LUT 28 so that for each set of (R, G, B)„ values there is a 
corresponding set of (R',G',B',)„ values. Means to store the 55 
transform LUT 28 is included in the compiler 82. This LUT 
28 is the the same transform LUT shown in FIG. 1 discussed 
earlier. 

The correction factor generator 76 produces the dR, dG 
and dB correction factors as follows. The input to the 60 
correction factor generator 76 is the difference signal pro- 
duced by the comparator 66, specifically [(Lab)„-(L*a'b')J. 
The individual components of this signal are: dL which 
equals (L — L*), da which equals (a — a') and db which equals 
(b— b*). The correction factors dR, dG and dB, and the 65 
difference signals dL, da, and db are related by the following 
relationships: 



(dbtdG)da - OafbOdb = 0) 
lOfr/ao ♦ (dafBR) - &a/BO * idbfBR)]dR + 

(Ob/aC) * ida/dB) - {da/BG) * mBB)\dB 



(BafdOdL - QUBC)da = 



(2) 



lida/BG) * (BUBR) - (BUBC) * (Ba/BR)\dR i 



|(9fl/3G) ♦ {BUdB) - {BUBO * {daJBB)\dB 
dL = (dLfBG)dG + {BUBR)dR + (BUdB)dB. (3) 

The actual numeric values for 3L/9R, 3L/3G, 3L/3B, etc. 
for use in solving the above system of equations, are 
obtained by the process schematically illustrated in FIG. 5. 
The input R„, G„,B„ values are changed by holding two the 
same and incrementing the third by one unit in incrementing 
module 96. This module comprises a memory for tempo- 
rarily holding the input values of R„, G„, and B„ and an 
arithmetic means for incrementing each of those values by 
1 and for applying various combinations of R„, G„, B„, 
R„-i-l, G„-i-l, and B„+l sets of values to the input of model 
42. 

An example which uses illustrative RGB values, will be 
used to explain the derivation of the numerical values for the 
partial derivatives needed to solve the above equations, (1) 
(2) and (3). Let the input R„, G„, and B„ values for a set of 
(R, G, B)„ values be 100, 100, 100. Assume that when these 
values are applied to both models 40 and 42, two sets of Lab 
values, (Lab)„ and (L'a'b')„ are produced such that the error 
E as previously defined is different than "0" or a preselected 
acceptable minimum. Incrementing module 96 produces 
over Hne 65 three new color value sets, R„, G„, (Bn+1), 
having the numerical values (100, 100, 101) in this illustra- 
tive example, R„, (Gn+l),Bn, having the numerical values 
(100,101,100) and (Rn+l), Gn, Bn, having the numerical 
values (101,100,100). These new values are applied to 
model 42 over line 63'. The resulting variations in the 
components of the (L'al}') n values produced by model 42 
are directed to an arithmetic calculator 82 over line 75 and 
to buffer memory 84. The calculator 82 performs the opera- 
tions: (I^L')/R-CR+1). (L— L')/G-(G+1), (L~L')/B-(B-h 
1), (a— a')/R-(R-hl), (a— a')/G-(G-Hl), (a— a')/B-(B-^l ), 
(b— b')/R-(R+l), b— b7G-(G+l) and b— b'/B-(B+l), to 
derive the partial derivative values (3L/3R, 3L/9G, etc.) used 
in equations (1), (2) and (3) above. 

The results of the calcidator 82 are supplied over line 33 
to another calculating unit 86. Tlie dL, da, and db values 
from comparator 66 are also directed to this calculating unit 
86. Unit 86 solves the system of equations (1), (2) and (3 ) 
for dR, dB and dG to derive the correction factors. 

A problem is sometimes encountered in determining a set 
of values that will produce a best color match on two 
different displays for the same input set of color values. The 
problem is that a minimum error may be arrived at, which 
may not be an absolute minimum, but only what may be 
referred to as a local minimum, and there may be other 
combinations of values (R', G', B')„ which could produce a 
smaller even a zero error and a better color match. To test for 
this possibility, an optimization operation is performed. A 
buffer memory 90 holds the input (R, G, B)„ values, the 
corresponding (R', G', B') values and the error *'E„" resulting 
from the use of this combination. The next set of (R, G, 
B)„-Hi values is processed and it, the corresponding (R', G'. 
B')n+i values, and resulting error "E„" are also stored. The 
following (R, G, B) n-H2, (R', G', B')„^2 and "Eg" are stored, 
etc., until a set of (R, G, B)„+jc and (R'G'BO^^^ values for 
which the error "E '*=0 is encountered. Once this true "0" 



02/24/2004, EAST Version: 1.4.1 



5.579,031 



9 

error has been identified, a backtracking operation is per- 
formed to determine if a "0" error can be obtained by 
correcting the set of (R, G, B)„_^_i values using' as an input 
to the model 42 the ^, G, B)„„j^ set of values (which is the 
set for which a "0'* error was obtained) while model 40 5 
receives the (R, G, B)„_^_i set. This operation will produce 
a new set (Rl', GV, Br)„_^i of values and a new error 
"El"„^.» for the input values (R, G, B)„^_i. If "El"„^ 
i="0" the (Rl', Gl*, Bl'),p^i set of values is substituted for 
the (R'GBOn set of values, and the transform LUT 28 is lo 
recompiled using these (Rl', Gl', Bl')^_^_^ values as corre- 
sponding to the (R. G. B)„^j set. 

If no **0" error is obtained with the above process, then the 
eariier arrived at result minimum is compared and if lower 
than the newly found minimum, the result slays Otherwise 15 
the R', G', B' values giving the least error are used in 
compiling the transform LUT 28 as they will result in a more 
accurate match. The optimization process is then repeated 
for the next (moving backwards) set of (R, G, B)^^^^ i^P^^ 
values and sequentially for all (R, G, B)„ sets of values for 20 
n<x always using the values from the set with *'E"=0 whose 
ny number is the closest to the tested set n number. When all 
stored sets of (R, G, B)„ values have been tested in this 
maruier, the memory contents of the buffer memory 90 are 
erased and a new set of input (R, G, B)„ beginning with the 25 
next n, i.e., n=x+i is applied over line 54 and the process of 
compiling the transform LUT 28 continues. 

For illustration purposes, lets assume that an input set of 
values (R, G, B)2o =(100. 100. 100) is applied to both the 
first model 40 and the second model 42 and that this results 30 
in a corresponding set of values (R',G',B')2o=(100. 102, 
103,) which gives a minimum acceptable error but different 
from zero. Let us further assume that the next set of input 
values is (R. G, 6)31 and that (R. G, 6)31 =(110. 100, 100). 
Also let us assume that the corresponding set of values (R', 35 
G', B')2, are (105. 100, 101), and that the error in this 
instance is zero. 

When the zero error is detected, the optimization process 
begins, and no new input sets of (RGB)„ values are applied 
to line 54. 40 

The set preceding the set which had the zero error, i.e., set 
of input values (R. G, 6)^0 ==(100. 100, 100) is applied to the 
input of model 40. However, the set (R',G',B')2^ =(105, 100, 
101) is applied lo the input of the second model 42 (rather 
than the set (R, G, B)2o) and a new set of corresponding 45 
(Rl'Gl'Bl')2o values is arrived at. i.e., 101, 102. 103, which 
will either give a zero error or another minimum error. If the 
result is a zero error or if the new error is smaller than the 
previous, the new (Rl'Gl'Bl')2o values are substituted for 
the original values (R'G'B')2o that were used earlier in the 50 
compilation of transform LUT 28 and will produce a closest 
match. 

The process described is best performed in a computer, or 
if in hardware, under the control of a computer suitably 
programmed to perform the above computations. The dia- 55 
gram of FIG. 5 is provided for illustration, rather than as an 
actual working embodiment of hardware sufficient to per- 
form the above process, so as to enable one to properly 
program a computer for this task. 

Referring again to FIG. 3, there is shown a block 70 60 
between the output of model 40 and the input to comparator 
66. This represents an optional color gamut correcting 
device, which may be used to match the output color gamuts 
of the two display devices. Gamut matching is often often 
used because the color range that is possible to reproduce 65 
with a particular display is usually not identical with the 
color range possible to reproduce with another. Thus, while 



10 

both displays have a large number of colors that are repro- 
ducible by both, there may be certain combinations of sets 
of color values that are outside the reproduction capability of 
one or the other display. While the use of the transform LUT 
28 in itself inherentJy provides color gamut matching 
between the two devices considered, additional color gamut 
matching may be used in the development of the LUT 28, by 
inserting a gamut matching device 70 as shown. This device 
will perform a preselected mathematical manipulation to the 
(Lab)n set of values so that there will not be any sets for 
which there are no corresponding (L'a'b')„ values from the 
model 42. Such mathematical manipulation may be as 
simple as simple truncation or as complicated as data 
compression. 

These types of mathematical operations are well known in 
the art of color matching. For a brief discussion of gamut 
matching and the necessity of color compression, also see 
• COLORIMETRY FOR ELECTRONIC IMAGING SYS- 
TEMS. An intensive short course for Scientists and Engi- 
neers", Copyright 1991, Munsell color science laboratory, 
Rochester Instimte of Tbchnology, by R. S. Bems and R. S. 
Hunter. Chapters 27 and 28. 

FIGS. 6, 7 and 8 illustrate the case where the images that 
arc to be matched are images produced in color proofcrs or 
printers which use a typical four color reproduction system 
with four color inks, such as Yellow (Y), Magenta (M). Cyan 
(C) and Black (K). 

Referring to FIG. 6, as in the previous example, a work 
station 10 generates digital color values representing pixels 
of a color image. Since a four color system is used and the 
reproducing equipment accepts YMCK signals, the work 
station 10 output is a set of YMCK signals over output line 
116. The signals are diTecte4 to a first printer 112 over line 
118 and to a second printer 114 over line 120. Ahead of the 
input to printer 114 is placed an adaptor 122 which again 
comprises a transform LUT 128, a receiving and determin- 
ing element 126 and an interpolator 130. Hie LUT 128 is a 
transform LUT which correlates a plurality of YMCK values 
to a second plurality of Y'M'C'K' values such that colored 
pixels produced using a set of YMCK values in printer 112 
and a corresponding set of Y'M'C'K' values from the LUT 
128 in printer 114 will be identical or virtually identical to 
the standard CIE average observer when viewed under 
similar illumination and surroundings. 

As for the previous example the interpolator 130 is used 
to supply corresponding values for YMCK input values not 
listed in the transform LUT 128. The transform LUT 128 has 
been compiled using substantially the same process as for 
the transform LUT 28 used in the previous example. As 
illustrated in FIG. 7, the work station 10 may generate a 
plurality of digital YMCK values such as shown in the 
following Table II: 

TABLE n 

Y: 0, 13, 26, 51, 76, 102. 128. 153. 178, 204, 230, 255. 
M: 0, 13, 26, 51, 76, 102, 128, 153, 178, 204, 230, 255. 
C: 0, 13. 26, 51, 76, 102, 128. 153. 178, 204, 230, 255. 
K: 0, 13, 26. 51, 89, 128, 178, 255. 



Each combinarion of Y,M,C,K. values from this Table 11 
represents a set (Y,M,C,K) n of color values. In this illus- 
u^tion, there are 13824 possible such sets of Y,M,CJC, and 
n=l to 13824 values. 

The 13,824 patches resulting from the different color 
value combinations above are printed by both printers 112 
and 114 to provide two sets 132 and 134 of 13,824 patches 
each. Again as before, a colorimetric measuring device 36 is 



02/24/2004, EAST Version: 1.4.1 



5,579,031 



11 



12 



used to read all of the above patches 132 and 134 and the 
colorimeter output is used to compile two LU'R representing 
the transfer fiinaions of the two printers 112 and 114. Again 
it is preferred that the colorimeter output is in the Lab color 
space and the LUTs correlate (Y.M,C,K)„ inputs to (L,a,b)„ 
outputs for each of the printers 112 and 114. These LUTs are 
used to construct models 140 and 142 using the generic 
model diagram shown in FIG. 4 for the two printers 112 and 
114, respectively. 

FIG. 8 illustrates the process followed for the compilation 
of the transform LUT 128. This process is similar to the 
process described earlier for the compilation of the trans- 
form LUT. The LUT, however, correlates two pluralities of 
color sets, each set having four rather than three variables. 
Of the four C„, Y„, M„,K„ variables, the three C„, Y„,M„ 
variables are independent variables, but the fourth, K„, is 
not, K„ is dependent on the other three from which it is 
derived. The relationship between CYM and K is well 
known in the graphic arts. 

Preselected sets of (C,Y,M,K)^ values from the work 20 
station 10 or elsewhere, are directed to the inputs of both 
model 140 and 142 over lines 154, 162, 156 and 158. It is 
preferred that the preselected values be the same as the 
values used to compile the LUTs used in the models 140 and 
142 and for purposes of this description the same value 25 
combinations resulting using the values in Table II. 

Each input set of values (C,Y.M,K)„ is applied directly to 
the input of model 140 over lines 154 and 158. The same set 
of values is also applied over line 156 to one input of adder 
160. The adder output is directed to switching circuit 198. 30 
The output of switching circuit 198 is controlled by a signal 
over line 194 and may be directed to the input of model 142 
or to transform LUT compiler 182 over lines 162 and 180, 
respectively. The input set of values is also applied to 
compiler 182 over line 155. 35 

The adder 160 has a second input on which are applied 
correction factors dC„, dY„ and dM„ produced in the cor- 
rection factor generator 176. When the (Y.M,C,K)„ set of 
values first appears at the adder 160 input, there are no 
correction factors applied to the adder and the values Y„,M„, 40 
C„.K„ are applied to the input of model 142 unaltered. 
Model 140 produces an output set of values (L,a,b)„ and 
model 142 produces another set (L',a',b')„. These two sets are 
compared in comparator 166. For this description, it is 
assumed that an optional gamut adjustment module 170 is 45 
not used to modify the (L,a,b)„ output of model 140 prior to 
applying it to the comparator 166 over line 172. 

The comparator 166 produces a difference signal [(Lab)„- 
(L'a*b')J. Tliis output is directed over line 174 to correction 
factor generator 176. Using this difference signal, correction 50 
factor generator 176 produces three correction factors dY„ 
(1), dM„ (ij, dC„fi) (in a manner similar to the manner in 
which correction factor generator 76 produced the dR. dG, 
and dB factors as previously explained). The correction 
factors are sent over line 178 and added to the corresponding 55 
components of the CY,M.C,K)„ set of values in adder 60, 
resulting in a new set of color values Y'„(i3=Y„+dY„(tj, 
M;(,j=M,-KiM„(i). and C^^,^C„-^C„^,^ on line 162. 

The fourth value, K„ is stored in element 177 which 
comprises a buffer memory and a means to modify the stored 60 
value K„ by predetermined increments. At this time K„ 
remains constant and these new Y'„(,), M'^^j, C'^^jj, and K„ 
values are again applied to the input of model 142 which 
produces a new output set of l^^2\* values on line 164 
which is again compared in comparator 166 with the set of 65 
(Lab)„ values from model 140. A new set of correction 
factors dY„(2). dM„(2)» and dC„(2) is produced and added to 



Y'nCD' ^nci) ^nd) produce a new set of Y„', M„', and 
C„' values. Again, an error "E* is defined as: 

E=((L-Lf-Ka-a')^+{b-b')2|'« 

The correction factor generator 176 also includes means 
to calculate, compare and store this error E„ and the corre- 
sponding (YMCK)n and (Y'M'CK)„ set of values that 
resulted in this error E„. 

The above process is repeated continuously correcting the 
original input values Y„,M„, C„, with K„ held constant until 
the error E is zero or a predetermined acceptable minimum. 
If a zero error is found, then the two sets of values produce 
a perfect match and the search terminates. Switch 198 is 
actuated and the two sets of values are used in the compiler 
182 to compile the transform LUT 128. If only a minimum 
is arrived at it is again stored in the buffer memory and K„ 
is incremented in unit 177 and the process repeated in search 
of another set of (Y"M"C"K')n values yielding another 
minimum. This new minimum is compared with the previ- 
ously found one. The (Y'M'C'K)„ values or the values 
(Y"M"C"BC')n which produced the best minimum are kept. 

The correction factor generator in this instance solves a 
system of differential equations that is of the same format as 
previously discussed. The equations only differ in terms of 
variable designation and are as follows: 

(db/dY)da-(da/dY)db= (la) 
[(dbtdY) * {da/dM) - ida/dY) * (db/dM)]dM + 

[(dbldY) * {dafdQ - OaBY) * idbfdC)\dC 

ida/BY)dL - (dUBY)da = (2a) 
[(dafdY) * (dUdM) - idUBY) • (da/dM)]dM + 

fOfl/an ♦ ouao - (dudo * {dafdouc 

dL = iWdY)dY + (BLfdM)dM + (dUdQdC 



(3 a) 

The partial derivatives are calculated in the same manner 
as the partial derivatives were calculated for the R, G, B, 
case, by incrementing the Y„,M„, C„ values by one, as 
explained earlier with the help of FIG, 5. 

The issue of whether a minimum error E„ is a true 
minimum error or a local minimum, meaning that there may 
be another smaller error resulting from a different combi- 
nation of (C, Y', M', K')„ values which could be found by 
starting with a different input set to the second model 142, 
rather than the (Y,M,C,K)„ set thai was used is again a 
problem. The same testing and optimization process 
described earlier can be used in this case too, to test and 
discover whether there may be other better combinations of 
(Y', M', C\ K')„ values, by holding in a buffer memory all the 
consecutive sets of (Y,M,C,K)„, the corresponding 
(Y*M'C'K)„ sets and error E„ as described earlier. However, 
it has been found that more accurate results may be obtained 
by following a somewhat more calculation intensive pro- 
cess. This process utilizes two input sets of values which 
have resulted in a zero error as inputs to obtain two opti- 
mized sets of (Yl', Ml', CI', Kl')„ and (Yr, M2', C2\ K2') 
values in addition to the (Y', M', C\ K\ set and then 
selecting the set that produces the best match, i.e., least error, 
from the three sets. 

Because of the presence of the fourth color K, which is not 
an independent color as the Y,M,C, and which by its 
dependency on the other three can introduce a number of 
color combinations resulting in matching colors, the exist- 
ence of two different sets of values both producing a perfect 
match is also possible. Further, the optimization process 
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must be able to rationally select the set that will produce the 
better overall match if used in the compilation of the 
transform LUT 128. When two sets produce perfect 
matches, then the one that produces better interpolation 
results for input values that are not present in the LUT will 
be selected. 

The optimization process adopted for this case, is based 
on selecting a predetermined fixed number "x" of sets of 
values (Y;M,C,K)„, the corresponding (V, C\ M\ K\ and 
the associated enor and, then, proceeding with the 
optinuzation calculation using zero error resulting set values 
from either side of a target set to derive two sets of possibly 
lesser error resulting values in addition to the set originally 
found. Thus, the input set (Y,M,C,K)„8, for instance, for 
which the set of (Y'C'M'K')„g gave us originally a minimum 
but non zero error, will be optimized using preceding and 
following values, to obtain (Yl', C\\ Ml', Kl')„g values with 
data from a set (V, C. M\ K')„y and (r2,M'2,C'2,K'2)„8 
using data from a set (Y' ,C',M',K')^ where y>8 and 2<8 and 
where both the "ny" and "nz" sets are the sets which 
produced a zero error closest to the set n-8. 

The following illustrative example is provided to explain 
the optinuzation process based on a fixed number of value 
sets discussed above. 

Assume that the transform table 128 is to be constructed 
using the (YMCK)„ sets of color values from Table II. Let 
us further assume that we are at a stage in the process where 
the sets of (YMCK)„ shown in Table EI have been sequen- 
tially introduced over line 154 to model 140, line 156 to 
adder 160 and line 155 to LUT compiler 128. As a result, 
corresponding sets of (Y'M'C'K') of values and associated 
error En values have been obtained. The errors En at this 
time have not been tested for optimum error condition, for 
enrors different than zero. The number of sets used for 
optimization will be limited to four (that is x=4), to make the 
illustration simpler. The four sets of values and associated 
error values held in memory at the beginning of the opti- 
mization are shown in Table HI. 



TABLE m 



Scl# 


(Y, M, C, K) 


(V, M'. C\ K) 


E 


n 


76, 76. 76. 89 


73. 72, 78. 90 


0 


n + 1 


76, 76, 102. 89 


73, 73, 100, 89 


0.1 


n + 2 


76. 76, 128, 89 


73. 75, 131, 91 


0 


n + 3 


76. 76. 153, 89 


73. 71. 152, 93 


0 



(Y*,M'.C in adder 160. As before, a new set of L'a'b' 
values is produced and again compared in the comparator 
166, and the process repeated, including the K incrementing 
step, until a minimum error value is found again. Assume 
these new values to be Y'l„^i=73, M*l^i-74, Cl„+i=lOO 
and K'l„+i=88. Let us further assume that this set of values 
results in a perfect match and a zero error. We temporarily 
store these results. 

Next we repeat the above process using this time the set 
of (Y'M'CK')„^.2 values as an input to the adder 160, while 
the (Y,M,C,K)„+i set continues to be applied to model 140. 
This time we find that another set of (Y*2, M'2, C'2, K7)^j 
values, i.e.. Y7^i=71, M'2„^i=76, C2„^t=100 and K'2„^,= 
91, where the error is again zero. These values are also 
stored. These result are shown in Table IV: 
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For input sets n, n+2 and n4-3 the error is zero. Thus, we 
have perfect match and can accept the (Y,M,C,K) and 
corresponding (Y", M', C, K') values in the compilation of 
the transfer table 128. However, we need to check the results 50 
of set n+1 to see if a better match may be available. 

We proceed as we did earlier in the optimization process 
discussed for the (R. G, B) example. However, the process 
is done with two sets of input values to model 142, the sets 
that gave a zero match located on either side of the set that 55 
is optimized. In this example, the set that is optimized is set 
n4-l and the two sets, that gave zero error on the two sides 
of set n+1, are the n set and the n+2 set. We apply the 
(Y,M,C.K)„+i (i.e., 76, 76, 102, 89) set of values again to 
model 140, but we apply to model 142 through die adder 160 60 
first the (Y\M',C',K')„ (i.e., 73, 72, 78, 90) set of values 
which is the set corresponding to the * V input, the set which 
gave a zero error in matching the input values for "n". 

The (Lab)„+j and (L'a'b')„ outputs from the two models 
are compared in comparator 166 and new correction factors 65 
dY'„, dM'„, dC'„, are generated which result in a new set, 
Yl„+i, Ml„^.i, C'„+,. Kl^i of values when added to the 



TABLE IV 


Sct# 


(Y. M, C. K) 


(Y', M', C, K') 


£ 


n 


76. 76. 76. 89 


73, 72, 78, 90 


0 


n+1 


76. 76, 102, 89 


(a) 73. 73, 100, 89 


0.1 






(b) 73, 74, 100, 88 


0 






Cc)71, 76, 100. 91 


0 


n + 2 


76, 76, 128. 89 


73, 75, 131. 91 


0 


n+3 


76, 76, 153, 89 


73. 71. 152. 93 


0 



The question now becomes which set of values (Y', M', C", 
K')„+i should we keep in compiling our table 128 for the 
input set (Y,M,C,K)„^.]. Obviously set (a) is the worst match, 
and we can discard it. Of the other two, it is prefened to keep 
the set that will provide the most accurate interpolation 
results, when interpolation involving those values must be 
used. One way to do this is to select the set that has values 
best distributed between the preceding "n" set and the 
following "n+2" set. A way to make this selection is as 
follows. Let Table V represent the two possibilities of 
selection. 

TABLE V 
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Sel 


Possibility U\ 


Possibility #2 


Input 


# 


CY', M', C. K) 


(Y\ M', C. K') 


(Y, M, C, K) 


a 


73, 72, 78, 90 


73, 72, 78, 90 


76, 76, 76, 89 


n + 1 (b) 


73. 74, 100, 90 




76, 76, 102, 89 


n + 1 (c) 




71. 76, 100, 88 


76, 76, 102, 89 


n + 2 


73, 75/131,91 


73, 75, 131,91 


76, 76, 128. 89 



We next compare the point spread deviation between 
input (Y,M,C,K) sets and the corresponding (Y',M',C'4C) 
sets. For instance, comparing the inputs between sets n and 
n+1, we have: Y.,i-Y„=0. M„,,-M„=0, C^rC„=26. K,,- 
K„=0 and comparing the inputs between sets n+1 and s+2, 
we have: Y^2-Y^,=0. M,,j-M,^,=0, C„^rC^i=26. K„.2- 
K„^.]=0. Next, we do the same using the values for Y', M', 
C, K' from the lint selection column (labelled "Possibility 
#1"). and we get: Y',^,-Y'„=0, M' =2, C'„^,-C =22, K'„«- 
K*„=-l . Thus, we see that for a zero change in the Y input 
values, there is a zero change in the Y' values. For zero 
change in the M input values, we have a change of +2 points 
in the M' value. For a 26 point change in the C value input, 
the C value changes by 22 points, or a short fall of -4. And 
for a zero input change in the K value we have a -1 point 
change in the K'. 

Similarly the differences observed between the n+2 and 
n+1 sets, yields a difference for Y=0 points, M^l, 0^5 
and K=+2. 

We will refer to these differences in the values as the 
deviation from the ideal situation where the increments of 
change between input values are the same as the increments 
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of change between corresponding values, that is where again compiled using the process illustrated in FIG. 8. The 

(C„-C„_i)=(C',,-C„_i). The average of the deviation magni- transform LUT 328 converts input YMCK values to output 

tudc AD„ for the deviation between sets n and n+l for Y', M', Y'M'CK' which when sent to converter 325 generate RGB 

C, K' from T^ble V selection #1 column is ADi=7/4 or 1 .75; values such that when displayed on the CRT display 314 will 

for sets n+1 and n+2. again for selection #1 the average is 5 produce an image having the same or substantially the same 

AD2=2. color appearance as the image produced by the color proofer 

For values from column selection #2 in Table V, the same 312 using the original YMCK values directly, 
process yields AD3=2.75 for sets n and iHl, and 2.25 for In discussing this invention there has been constant ref- 
AD4 for sets n+1 and n+2. Thus, after comparing the average erence to "interpolation". Interpolation leading to derived 
deviations for the two selections, set n+l(b) in selection #1 10 output sets of color image values based on input sets of 
is seen to yield set values that have magnitudes which, when image values is used in performing the process of this 
taken as adjacent pairs of values with the preceding or invention. It is evident that if interpolation is not used, the 
following sets of values define color spans that on the LUTs employed both in the models and in the adaptors must 
average deviate less from the spans defined by the corre- supply values for all possible combinations of input color 
sponding input color pairs than the values obtained from the 15 values. In a digital system operating in an 8-bil domain, for 
sets of values from selection #2. Therefore, the set of (Y', M', a three variable system, i.e., RGB, the possible input value 
C, K*)^i(j,j values is used in the transform table compila- combinations are 256*256*256 or 16,777,216 combina- 
tion. If set (a) could not be discarded as obviously the worst tions. To develop the necessary LUTs, each of these corn- 
match, it could be compared to sets (b) and (c) like set (b) binations would have to be printed as a patch, in accordance 
was compared to set (c) above to determine which set is the 20 with this process, measured, compared, and optimized, 
best match. clearly an impractical process because of the shear magni - 

FIG. 9 shows a situation where image color matching is tude of numbers and calculations involved. Thus, limited 

desired in images produced with a color printer 214 which numbers of sets are selected as illustrated in Tables I and II 

accepts RGB values (through a converter 225 that converts and interpolation is used for intermediate sets of input color 

RGB inputs to YMCA using a mathematical relationship that 25 values. 

the primer manufacturer has determined as an accurate color The particular interpolation process used may vary, 
conversion) and a CRT monitor 212. For this arrangement, Depending on the accuracy one wishes in the system more 
as shown in FIG. 9, the output of work station 10 is an RGB or less calculation intensive interpolations may be used. In 
signal over line 216. It is used to drive the CRT display 212 practicing this invention, it is preferred that a less calculation 
directly and the color printer 214 through a 3 to 4 color 30 intensive interpolation method is used in the model inter- 
conversion unit 225 which converts RGB inputs to YMCK polators, while a more sophisticated method is used in the 
outputs; an adaptor 222 comprising a transform LUT 228 adaptor interpolators. Examples of both types are given 
constructed in accordance with this invention is placed below. However, any of the interpolation methods described 
ahead of conversion unit 225. herein below, as well as other interpolation formulae, such 
This LUT 228 is compiled using a process similar to the 35 as quadratic interpolation, can be used for either the adaptor 
process illustrated in FIG. 8. However, in generating the interpolators, or the model interpolators, the final choice 
model 44 for the printer 214, the converter unit 225 is being based on desired accuracy of the matching process, 
included with the color printer 214 as an integral part and the availability of resources. This invention is not 
thereof. Thus, this model 44 which is schematically shown restricted to the use of one interpolator for adaptors 22, 122, 
in generic format in FIG. 11 has an input 45' which leads to 40 222, 322 and a different one for the models 40, 42, 140 and 
converter 47. The output of converter 47 over line 45 is 142 even though such arrangement is both used and pre- 
directed to the input of a receiving and determining element ferred. 

46 and therefrom in a similar arrangement to that shown in The preferred interpolation performed in interpolators 30, 

the generic model of FIG. 4 produces an output 52 either 130, 230 and 330 shown in the adaptors shown in FIGS. 1, 

direcUy through LUT 48, or by interpolation through inter- 45 6, 9 and 10, respectively, tiiat is the interpolator 30, 130, 230 

polator 50. or 330 that accompanies a transform LUT 28, 128, 228 or 

In the arrangement shown in FIG. 10, the output of the 328, respectively, is tetrahedral interpolation performed 

work station 10 is a YMCK signal which drives directly through a programmed computer. A pair of appropriate 

color proofer 312 over lines 316 and 318. The YMCK output programs for programming such a computer to perform this 

is also used to display an image on a CRT monitor 314. This 50 interpolation for three and four variables are included with 

monitor 314 is able to accept YMCK input values because this application as Appendices "A" and "B" which appear 

of a 4 to 3 color conversion that occurs in a color converter inmiediately before the claims, 

325 prior to being used by the monitor 314. To obtain more Linear interpolation, tri-linear or quadralinear depending 

accurate image color matching, an adaptor 322 which on the particular color space in which the color values exist, 

includes a receiving and determining element 326 and a 55 Is the preferred choice of interpolation in models 40, 42. 140 

transform LUT 328 for converting YMCK input values to and 142 shown in FIGS. 3 and 8, used in deriving the 

Y'M'CK' is placed ahead of the converter 325. As before an transform LUTs 28. 128, 228. 328, as it is simpler, and often 

interpolator 330 is included in the adaptor to derive values adequate to produce good results. The algorithm for this type 

not present in the LUT 328. The Y'M'CK' values are of interpolation as applied in this invention is next explained 

directed to a CRT display 314 over line 324 through existing 60 for the more complex case of a four to three transform, i.e., 

converter 325. from (YMCK) input values to (Lab) as would be used in 

The transform LUT 328 is also generated using the models 140 or 142. The three to three (RGB color space) and 

process described earlier, wherein two models 44 are created four to four (YMCK color space) transform cases are 

representing the two display devices 312 and 314. Hie CRT similar. 

display device model 314, however, as in the previous 65 The numerical values for Y,M,C,K given in Table U 

example, includes the 4 to 3 convener 325 and is also similar represent data points along a system of axes Y, M, C, and K. 

to the generic model illustrated in FIG. 11. The LUT 328 is In the case of three variables, i.e., R.G,B color values, this 
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may be an orthogonal axis system that may be easily 
visualized. With four variables visualization is very difficult, 
however, such actual visualization is not necessary to per- 
form the interpolation. 

The data points representing predetermined color values 
of a variable along the axis for this variable are separated by 
a specified number of units representing a distance measured 
along this axis. Thus, the data points Y=102 and Y=128 are 
separated by a distance of 26 units along the Y axis, while 
the data poinls M=128 and M=153 are separated by 25 units 
along the M axis. An input value may fall either on a data 
point, or between data points. When it falls on a data point 
the LUT provides a corresponding value. If the input value 
fall between data points, then an output is calculated as 
follows. 

Let be the ratio of the distance of the input value from 
a preceding data point along an axis to the distance between 
that data point and the next consecutive data point along this 
axis. Thus, if the color value set has four components, (i=0. 
1, 2, or 3, for the Y, M, C, K components) and Y„ is a data 
point along the Y axis, Yn+l is the next data point along this 
axis, and Y, is the Y component value of a color input set 
of values whose numerical value lies between Y„ and Y„^.i, 
then: F.-(Y1-Y„)/(Y„,,-YJ for i^. 

We also define a quantity Qi=l-Fi, and the coefficients tj 
(where j=l, ... 16) as follows: 



ti^Qo*Fi^F2*Fi 
ti = FQ*Qi*F2*Fi 
U = FQ*FrQ2*Fi 

/l6 = G0*Cl*fi2*C3 



t6 = QQ''Ql*F2*Fi 

t7-'QQ^Fx*Q2*Fi 
t9 = Fo*Qi*Q2*Fi 



hx=^Fo*Fi*Qi*<h 
t^l = Q^*F^*Ql*Qi 

lH = QQ*Ql*F2*<h 



10 



18 

'3 * ^nM-/,*fti,ow-j,j&M-/j + '4 * lY»i+;,AftH-;,oi,»>+/; + 

'11 * L(Y^jj^f^^.i^Cm.Xit) + '12 * £'(]WJ,Wn.Oi.jen> + 
/1 3 * iflXAfw+i.OUI&iJ + '14 * L(YH^UH.Cm^J,Kn) + 

Similarly, 

'11 * 0(Yn+J,Mn¥t.Qi.Kn} + '12 * Of C«,ICit J + 

'13 * a(m,Mn+i,Ctt.lCn) + '14 ♦ a(YH,Mn.at+/.Kn} + 

'15 * 0(yj^f^cn.Kni^l) + '16 * Ofrn.Mfi.Ch.Kr)' 

And, 

30 bj=ti* bfYn+t,Mn+I,Cn+J,Kn^i) + '2 * ff(Yt,.MtnJ.Cm^t.JCn*I) + 

'3 * b(Yn^j^Mn,Cni^l.Kn*Jl + '4 * b(y„^tj^^t,Ch,Ki>*Jl + 

'5 * i'frB+;,Mn+i.OH-/.»ij + '6 * *r>XM«.ai+/.A«w> + 
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Let l^(^Yn,Mn,CnJCn)^ ^YnMn^CnJCn^* ^(yn,Mn,Cn,Kn)* repICSCnt 

the Lab values from the LUT for the model 140 for instance, 
for an input Y„.M„.C„.K„, Uyr^iMn.cn^Kn)* a (K„+i.Mn.oi.A:J. 
\yn^iMn,cn.Kn) ^ab valucs for an input Y„+i, M„, C„, and 

^/i' ^irnMnfl,CnJCny ^Yn,Mn-^l,Cn,Kn)* ^YnMn+UCnJCny 

values from the LUT for an input Y„, M„^.i, C„ and K„, and 
so on for all 16 combinations of input values (Y, M, C, K)„ 
and (Y, M, C, K)^i. The Lab values Lj, a„ and bj for an 
input set of Y1,M1,C1,K1 color values falling between the 
n and n+1 sets of values in the LUT, arc given by the 
following equations: 



40 



45 



'9 * b(YK^JJ^^„,Cm.Jb*^■J) + flO * b(Yn^riMn.Cn*t.Kn) + 
'11 * b(ym+l.M>t*i.Ch,Kt>} + '12 * fc^y>,+/.Jlto.CW.J&i> + 
'13 * b(Yn.Mn-*^J.CKKnt + '14 * iffl^JWrvOi+AKh) + 

'15 * fcfr«,AfivCHJCm-/> + '16 * bfYr.Mn.Q,.Xmy 

Those skilled in the art, having the benefit of the teachings 
of the present invention as herein above set forth, can effect 
nimierous modifications thereto. These modifications are to 
be construed as being encompassed within the scope of the 
present invention as set forth in the appended claims. 
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APPENDIX A 



1 /*******«****«******«*****♦****«************************************* 

2 

3 (c) Copyright 1992 

4 E.I. DuPont de Nemours and Company { Inc , ) 

5 Wilmington^ Delaware 19898. 
6 

7 Module: tr3us_t.c 
8 

9 Purpose: 
10 
11 
12 
13 
14 
15 

16 Author: C. B. Chittineni 

17 

18 Date Create: June 1, 1992 
19 
20 



21 Usage: tr3us_t {n,ns,xl,x2, x3/ 1^ x, y, z, val) 
22 

23 where n - no of outputs (int short) 

24 ns - total number of data sample points per output (int short) 

25 xl - data point with in the cube along xl axis (int short) 

26 x2 - data point with in the cube along x2 axis (int short) 

27 x3 - data point with in the cube along x3 axis (int short) 

28 1 - length of cube (int short) 

2 9 X - data vector at current node (least rapidly varying) 

30 (unsigned char) 

31 y - data vector at current node (second least rapidly 

32 varying) (unsigned char) 

33 2 - data vector at current node (most rapidly varying) 

34 (unsigned char) 

35 val- pointer to interpolated values (n) (int short) 
36 

37 Modification: (Chittineni: July 9, 1992)Changed the vertex numbering 

38 as (n) -> (X y z) (0 - 1 1 1, 1 - 0 1 1, 2 - 1 0 1, 

39 3-00 1, 4-11 0, 5-01 0, 6-10 0, 

40 7-000) (Liang's data files format) . 
41 

42 / 
43 

44 tr3us_t(n, ns, xl, x2, x3,. 1, x, y, z, new) 

45 short new[3], n, ns, xl, x2, x3, 1; 
4 6 unsigned char *x, *y, *2; 

47 { 
48 

4 9 int short wO, wl, w2, w3, a, b, c, d; 

50 int itm, j, jl, dt; 

51 int j2, j3, j4; 

52 unsigned char *xO, *yO, *zO; 
53 

54 itm=0; 

55 if (x2 >- xl) { . 

56 if (x3 >= x2) { 

57 itm=3; 

58 lelse if (x3 >== xl) { 

59 itm=5; 

60 )else{ 

61 itm=^4; 

62 } 

63 )else if (x2>-x3) { 



A function for tetrahedral type highspeed interpolation 
for performing 3 to 3 or 3 to n color coordinate 
conversions for the case of uniform sampling along 
coordinate axes. The vertices of the cube are 
numbered sequentially in a natural binary order 
(x3 x2 xl). 
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64 itin=0; 

65 }else if (x3 >- xl) { 

66 itm-2; 

67 }else{ 

68 itin=l; 

69 } 
70 

71 switch (itm) { 

72 

73 case 0: 

74 wO=(l-xl); wl«(xl-x2); w2-(x2-x3); w3»x3; 

75 new[0]-(wO*x[7]+wl*x[3]+w2*x[l]+w3*xIO])/l; 
7 6 newtl]=(wO*y(7]+wl*y(3]+w2*y[l]+w3*ylO])/l; 

77 new[2]=(wO*2[7]+wl*z(3]+w2*z(l]+w3*ztO])/l; 

78 break; 
79 

80 case 1: 

81 wO=(l-xl); wl={xl-x3); w2=(x3-x2); w3=x2; 

82 newt0]=(w0*x[7]+wl*x[3]+w2*x[2]+w3*x[0]) /I; 

83 new[l]*(w0*y[7]+wl*y [3]+w2*y[2]+w3*y [0]) /I; 

84 new[2]'=(w0*z[7]+wl*z[3]+w2*z[2]+w3*z[0])/l; 

85 break; 
86 

87 case 2: 

88 w0=»(l-x3); wl=(x3-xl); w2=(xl-x2); w3"x2; 

89 new[0]-(w0*x[7]+wl*x(6]+w2*xl2]+w3*x[0]) /I; 

90 new[l]-(w0*y[7]+wl*y[6]+w2*y[2]+w3*y [0]) /I; 

91 new[2] = (w0*z[7]+wl*z(6]+w2*z[2]+>f3*2 [0]) /I; 

92 break; 
93 

94 case 3: 

95 w0=(l-x3); wl«{x3-x2); w2=(x2-xl); w3=xl; 

96 new[0]«(w0*x[7]+wl*x(6]+w2*x[4]+w3*x[0]) /I; 

97 new[l]-=<w0*y[7]+wl*y[6]+w2*y[4]+w3*yt0]) /I; 

98 new[2]-(wO*2[7]+wl*z(6]+w2*z[4]+w3*2[0]) /I; 

99 break; 
100 

101 case 4: 

102 wO=(l-x2); wl={x2-xl); w2=(xl-x3); w3«x3; 

103 new[01-(w0*x[7]-fwl*x[51+w2*xri]+w3*x[01)/l; 

104 new[l)=(w0*y[7]+wl*y[5]+w2*y[l]+w3*y[0])/l; 

105 new[2]«(w0*z[7]+wl*z[5]+w2*z[l]+w3*2[0])/l; 

106 break; 
107 

108 case 5: 

109 wO=(l-x2); wl»(x2-x3); w2=(x3-xl); w3=xl; 

110 new[0]=(w0*x[7]+wl*x(5]+w2*xl4]+w3*x[0])/l; 

111 new[l]-(w0*y[7]+wl*y[5]+w2*yt4]+w3*yt0])/l; 

112 new[2]=(w0*2[7]+wl*z [5]+w2*2 [4]+w3*2 [0] ) /I; 

113 break; 
114 

115 } 

116 return; 

117 } 
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APPENDIX B 

1 + it****** 

2 

3 <c) Copyright 1992 

4 E. I. DuPont de Nemours and Company (Inc.) 

5 Wilmington, Delaware 19898. 
6 

7 Module: tr4us_t,c 
8 

9 Purpose: A function for tetrahedral type highspeed interpolation 

10 for performing 4 to 4 or 4 to n color coordinate 

11 conversions for the case of uniform sait^ling along 

12 coordinate axes. The vertices of the hypercube are 

13 numbered sequentially in a binary order (xl x2 x3 x4). 

14 It is to be noted that the ordering is inportant 

15 in arranging the sampled data and in passing the 

16 arguments to the function. 
17 

18 Author: C, B. Chittineni 
19 

20 Date Create: June 1, 1992 
21 

22 Usage; tr4us_t (n,ns/Xl,x2, x3^ x4, l,cbl,rabl,ybl, kbl, val) 
23 

24 where n - no of outputs (int short) 

25 ns - total niimber of data sample points per output (int) 

2 6 xl - data point with in the hypercube along xl axis (int short) 

27 x2 - data point with in the hypercube along x2 axis (int short) 

28 x3 - data point with in the hypercube along x3 axis (int short) 

29 x4 - data point with in the hypercube along x4 axis (int short) 

30 1 - length of hypercube (int short) 

31 cbl- data vector at current node (least rapidly varying) 

32 (unsigned char) 

33 mbl- data vector at current node (second least rapidly varying) 

34 (unsigned char) 

35 ybl- data vector at current node (third least rapidly varying) 

36 (unsigned char) 

37 kbl- data vector at current node (most rapidly varying) 

38 (unsigned char) 

39 val- pointer to interpolated values (n) (int short) 
40 

41 Modification: (Chittineni : July 9, 1992) Changed the vertex numbering 

42 as (n) -> (cbl rabl ybl kbl) (0-llllr 1-0111, 

43 2-1011, 3-110 Ir 4-1110, 5-0011, 

44 6-010 1, 7-0110, 8-100 1, 9-1010, 

45 10- 1 1 0 0,11 - 1 0 0 0,12 - 0 1 0 0>13 - 0 0 1 0, 

46 14- 0 0 0 1,15 - 0 0 0 0) (Liang's data format) . 
47 

48 ******* ************* ************************************** *★★*★******/ 
49 

50 #include <stdio.h> 

51 finclude <math.h> 



52 tr4us_t(n, ns, xl, x2, x3, x4, 1, cbl, rabl, ybl, kbl, new) 

53 int short n, xl, x2, x3, x4, 1; 

54 int short *new; 

55 int ns; 

56 unsigned char *cbl, *mbl, *ybl, *kbl; 

57 { 
58 

59 int short wO,- wl, w2, w3, w4, a, b, c, d, e; 

60 int itm, j, jl, j2; 
61 

62 /* Locate the tetrahedran */ 

63 if (xl >= x2) { 
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64 if {x2 >• x3) 1 

65 if (x3 >- x4) { 

6 6 itn^l ; 

67 }else if (xl >«= x4) { 

68 if (x2 >= x4) { 

69 ititi=2; 

70 }else{ 

71 itm-3; 

72 } 

7 3 }else { 

74 itin=19; 

75 ) 

7 6 }eise if (x3 >» x4) { 

77 if<xl>=x4){ 

78 if(x2>-x4)< 

79 if(xl>=x3){ 

80 itm=5; 

81 lelse{ 

82 itm=13; 

83 } 

84 }else if (xl >= x3) { 

85 itrn=6; 

86 )else{ 

87 itm=14; 

88 } 

8 9 )else{ 

90 itin=15; 

91 ) 

92 }else if (xl >« x4) { 

93 itm='4; 

94 }else if (xl >= x3) I 

95 itm-20; 

96 }else { 

97 itm=21; 

98 } 

99 }else if (x2 >= x3) { 

100 if (x3 >" x4) { 

101 if{xl >= x4){ 

102 if (xl >= x3> { 

103 itin-7; 

104 }else [ 

105 itnt^ll; 

106 } 

107 }else{ 

108 ititt=12; 

109 } 

110 }else if(xl x4) { 

111 it:m=8; 

112 }else if (x2 >= x4) { 

113 if (xl >= x3) { 

114 itm=9; 

115 }else{ 

116 itm=10; 

117 } 

118 )else if(xl >= x3) { 

119 itm=23; 

120 )else{ 

121 itm=24; 

122 1 

123 }else if (x3 >- x4) { 

124 if (xl >= x4) { 

125 itm=17; 

126 }else if (x2 x4) { 
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127 itm=18/ 

128 }else{ 

129 itin=16; 

130 } 

131 }else{ 

132 itm=22; 

133 } 
134 

135 /* fprintf (stderr, "%d \n" , itm) ; */ 

13 6 

137 switch (itm) { 

138 

139 case 1: 

140 wO=(l-xl); wl-(xl-x2); w2=(x2-x3); w3-<x3-x4); w4-x4; 

141 new[0] = (w0*cbl(15]+wl*cbl(ll]+w2*cbl(10)+w3'^cbl[4]+w4*cbl(0])/l; 

142 new [1 ] = (wO*inbl [151 +wl*inbl [ 11 ] +w2*nibl [ 10 ] +w3*inbl [ 4 ] +w4*mbl [ 0 1) /I; 

143 newt2]=-(w0*ybl[15)+wl*ybllll]+w2*ybl[10]+w3*ybl[4]+w4+ybl[0])/l; 

144 new[3] = {w0*)cbl(15]+wl*kbl[ll]+w2*kbl[10]+w3*kbl[4]+w4*kbl[0])/l; 

145 break; 
146 

147 case 2: 

148 wO=-(l-xl); wl«(xl-x2)/ w2=(x2-x4); w3={x4-x3) ; w4=x3; 

14 9 new[a]='(w0*cbl[15)+wl*cbl[ll]+w2*cbl[10]+w3*cbl[3]+w4*cbl(0])/l; 

150 new[l] = (wO*inbl[15]+wl*nibltll]+w2*naDl [10]+w3*n±»l[3]+w4*inbl(0]) /I; 

151 new[2]-(w0*ybl[15]+wl*ybl[ll]+w2*ybl[10]+w3*ybl[3]+w4*yblI0])/l; 

152 new[3]-(w0*kbl(15]+wl*kbl[ll]+w2*kbl[10]+w3*kbl[3]+w4*kbll0])/l; 

153 break; 
154 

155 case 3: 

156 wO=(l-xl); wl=(xl-x4); w2=(x4-x2); w3=(x2-x3); w4=x3; 

157 new[0] = (w0*cbl[15]+wl*cbl[ll]+w2*cbl[8]+w3*cbl[3]+w4*cbl[0] )/l; 

158 new[l]«(w0*inbl[15]+wl*inbl[ll]+w2*iTa5l[8]+w3*iTt>l(3]+w4*mbl(03 )/l^^ 

15 9 newt2 3=(wO*ybl[15]+wl*ybl[lll+w2*ybl[81+w3*yblt3]+w4*ybl[0] )/l; 

160 new[33=(w0*kbl[15]+wl*kbl[ll]+w2*kbl[8]+w3*kbi[3]+w4*kbl[0] )/l; 

161 break; 
162 

163 case 4; 

164 wO-(l-xl); wl°(xl-x4); w2=(x4-x3); w3=(x3-x2); w4«x2; 

165 new[0]=(wO*cbl(a5]+wl*cbl[ll]+w2*cbl[8]+w3*cbl[23+w4*cbl[0] )/l; 

166 new[l]-(w0*inblll5]+wl*inbl[ll]+w2*rnbl[8]+w3*mbl (23+w4*inbl [0] ) /I; 

167 new[23«(wO*ybl|151+wl*ybl[113+w2*ybl[83+w3*ybl[2)+w4*ybl[03)/l; 

168 newt33=(w0*kblI15]+wl*kbl[113+w2*kbl[83+w3*kbl(23+w4*kbl[O3)/l; 

169 break; 
170 

171 case 5: 

172 w0=(l-xl); wl«(xl-x3); w2»(x3-x2); w3={x2-x4); w4=x4; 

173 - new[0)«(w0*cbl(153+wl*cbiIll]+W2*cbl[93+w3*cbir4]+w4^ " 

174 new[13==(wO*mbl[153+wl*inbl[ll]+w2*inbl[9lH^3*mbl(4]+w4*inbl[03 )/l; 

175 new[23 = (w0*ybl[15]+wl*ybl[ll]+w2*ybl[91+w3*ybl(4]+w4*ybl[03)/l; 

176 newE31-{wO*kbl(lS3+wl*kbl[ll]+w2*kbl[93+w3*kbl(43+w4*kbl[03)/l; 

177 break; 
178 

179 case 6: 

180 w0=(l-xl); wl=(xl-x3); w2=(x3-x4); w3«(x4-x2); w4=x2; 

181 new[03=(wO*cbl(153+wl*cbl[113+w2*cbl[93+w3*cbl(23+w4*cbl[01 )/l; 

182 new[13-(wO*inblll53+wl*ii±>l[113+w2*rnbl[9]+w3*rnbl[23+w4*mbl[03)/l; 

183 new[23 = {w0*ybl(15)+wl*ybl[113+w2*ybl[9]+w3*ybl(2)+w4*ybl[0))/l; 

184 new[33={w0*kblI15]+wl*kbl[ll]+w2*kbl[93+w3*kbl(2]+w4*kbl[0))/l; 

185 break; 
186 

187 case 7: 

188 w0=(l-x2); wl«(x2-xl>; w2-(xl-x3); w3-(x3-x4); w4-x4; 

189 new [03 = (w0*cbl [153 +wl*cbl [12] +w2*cbl [10] +w3*cbl [43 +w4*cbl [03 ) /I; 
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190 new[l]-(w0*inbl[15]+wl*n^l[12]+w2*mbl[10]+w3*na5l[4]+w4*mbl[0])/l; 

191 new[23-(w0*ybl[15]+wl*yblC12]+w2*ybl[10)+w3*ybl[4]+w4*ybi[0])/l; 

192 new(3]==(w0*kbl[15]+wl*kblll2]+w2*kbl[10)+w3*kbl[41+w4*kblt0])/l; 

193 break; 
194 

195 case 8: 

196 w0-(l-x2); wl-(x2-xl); w2-(xl-x4); w3-(x4-x3); w4-x3; 

197 new[0]-(w0*cbl[15]+wl*cbl[12]+w2*cbl[10]+w3*cbl[3]+w4*cbl[0])/l; 

198 new(l]-(w0*mbl[15]+wl*inbl[12]+w2*iribl[10]+w3*inbl[3)+w4*mbll0])/l; 

199 new(2]-<w0*ybl[15]+wl*ybl[12]+w2*ybl[10]+w3*ybl[3]+w4*ybl[0])/l; 

200 new(3]-<w0*kbl(15]+wl*kbl[12]+w2*kbl[10]+w3*kbl[3]+w4*kbll0])/l; 

201 break; 
202 

203 case 9: 

204 w0-{l-x2); wl-<x2-x4); w2=(x4-xl); w3-(xl-x3); w4-=x3; 

205 new[0]-(w0*cbl(X5]+wl*cbl[12]+w2*cblt6]+w3*cbl[3)+w4*cbl[0n/l; 

206 new[l)-(w0*na5l[15]+wl*itibl[12]+w2*na3l[6]4-w3*n0:)l[3]+w4*xt^^ 

207 new[2]=(w0*ybl[15]+wl*ybl{12]+w2*ybl[6]+w3*ybl[3]+w4*ybl[0])/l; 

208 new[3]=(w0*kbl[15]+wl*kbl[12]+w2*kbl[6]+w3*kbl[3]+w4*kblI0])/l; 

209 break; 
210 

211 case 10: 

212 wO-(l-x2); wl'^ (x2-x4) ; w2=(x4-x3); w3-(x3-xl); w4«xl; 

213 new[0]=(w0*cbl[15]+wl*cbl[12]+w2*cbl[6]+w3*cbl[l]+w4*cbl(01)/l; 

214 new [ 1 ] « ( wO *inbl [ 1 5 ] +wl *mbl [ 12 ] +w2 *inbl [ 5 ] +w3 *inbl [ 1 ] +w4 *nibl [ 0 ] ) / 1; 

215 newC2]=(w0*ybl[15]+wl*ybl[12]+w2*ybl[6]+w3*ybl[l]+w4*yblC0])/l; 

216 new[3]-(w0*kbl[15]+wl*kblC12]+w2*kbl[6]+w3*kbl[l]+w4*kbl[0])/l; 

217 break; 
218 

219 case 11: 

220 w0«(l-x2); wl«{x2-x3); w2=-(x3-xl); w3«(xl-x4); w4-x4; 

221 new[0]=(w0*cbl[15]+wl*cbl(12]+w2*cblt73+w3*cbl{4]+w4*cbl[0]) /I; 

222 new(l)-(w0*irtol[15]+wl*rttol[123+w2*mbl[7]4-w3*na3l[4]+w4*mbl[0])/ 

223 new[2]=«(w0*ybl[153+wl*ybl[12]+w2*ybl[7]+w3*ybl[4]+w4*ybl[0])/l; 

224 new[3]=(w0*kbl[15]+wl*kbl[12]+w2*kbl[7]+w3*kbl[43+w4*kbl[01)/l; 

225 break; 
226 

227 case 12: 

228 wO=(l-x2); wl={x2-x3); w2={x3-x4); w3*(x4-xl); w4-xl; 

229 new[0]=(w0*cbl[153+wl*cbl(12]+w2*cbl[7]+w3*cbl[l]+w4*cbll01)/l; 

230 new [ 1 ] - { wO *mbl [ 15 ] +wl*nibl ( 12 ] +w2*inbl [ 7 ] +w3 *nibl ( 1 ] +w4 *mbl [ 0 ] ) /I; 

231 newt2]-{w0*ybl[15]+wl*ybl[12]+w2*ybl[7]+w3*ybl(l]+w4*ybl[0])/l; 

232 new[3]-(w0*kbl[15]+wl*kblll2)+w2*kbl[7]+w3*kblll]+w4*kbl[0])/l; 

233 break; 
234 

235 case 13: 

236 w0=d-x3); wl=(x3-xl); w2=(xl-x2); w3=(x2-x4); w4-x4; 

237 new(0]={w0*cbl(15]+wl*cbl[13]+w2*cbl[9]+w3*cbl[4]+w4*cbl[0])/l; 

238 new [1] - (wO*robl [ 15]+wl*mbl [ 13] +w2*mbl (9] +w3*mbl [4] +w4*inbl [0] ) /I; 

239 new [ 2 ] = ( wO * ybl [ 1 5 ] +w 1 *ybl 1 1 3 ] +w2 *ybl ( 9 ] +w3 * ybl [ 4 ] +w4 * ybl 1 0 ] ) / 1 ; 

240 new(31-(w0*kbl[15)+wl*kbl[13]+w2*kbl[9]+w3*kbl[4]+w4*kbll0])/l; 

241 break; 
242 

243 case 14: 

24 4 w0-(l-x3); wl-(x3-xl); w2=(xl-x4); w3={x4-x2); w4=x2; 

245 new[01«(w0*cbl(15]+wl*cblll3]+w2*cbl[9]+w3*cbl[2]+w4*cbl[0])/l; 

24 6 new[l)-(wO*iabl(15]+wl*mbltl3]+w2*mbl[9]+w3'^nibl(2]+w4*mblt0]) /I; 

24 7 new[2]«(wO*ybl(15]+wl*ybl[13]+w2*ybll9]+w3*ybl[23+w4*ybl(01)/l; 

248 new[3]-(w0*kblll5]+wl*kbl[13]+w2*kbl(9]+w3*kbl[2]+w4*kbl[0))/l; 

24 9 break; 

250 

251 case 15: 

252 w0»(l-x3); wl=(x3-x4); w2=(x4-xl); w3=(xl-x2); w4=x2; 
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253 new[0]-(w0*cbl(15]+wl*cbltl3]+w2*cbl{5]+w3*cblE2]+w4*cbl(0]) /I; 

254 new[13 = (wO*^0^1(153+wl*^nbl{13]+w2*Inbl[5]+w3^^a^l[2]+w4^ 

255 new(2] = (w0*yblll5]+wl*yblll3]+w2*ybl[5]+w3*ybl[2]+w4*ybl[0])/l/ 

256 new[3]«(w0*]cbl[15]+wl*kbl[13]+w2*kbl[51+w3*kbl(2]+w4*kbl[0])/l; 

257 break; 
258 

259 case 16: 

2 60 wO=(l-x3); wl=(x3-x4)/ w2=(x4-x2); w3={x2-xl) ; w4=xl; 

2 61 new{0]«(w0*cbltl51+wl*cbi[13]+w2*cbl[5]+w3*cbl[l]+w4*cbl(0]) /I; 

2 62 new[13 = (w0*nibl[15]+wl*nibl[13]+w2*nibl(5]+w3*nibl[l)+w4*iTa^^ 

2 63 new[2]-<w0*ybl[15]+wl*ybltl3]+w2*ybl(5]+w3*ybl[l]+w4*ybl[0])/l; 

264 new[3] = (wO*kbl[151+wl*kbltl3]+w2*kblI53+w3*kbl[l]+w4*kbl(0]) /I; 

2 65 break; 

266 

2 67 case 17: 

268 wO-(l-x3); wl-(x3-x2); w2=(x2-xl); w3=<xl-x4); w4=x4; 

2 69 new[0]«(w0*cbl[15]+wl*cbi[13]+w2*cbl[7]+w3*cbl[4]+w4*cbl[0])/l; 

270 new[13-(w0*inbl[15]+wl*ria3l[13]+w2*mbl(7]+w3*riibl[4]+w4*rnb^ 

271 new[2]-(w0*ybl[15]+wl*ybl[13]+w2*ybl(71+w3*ybl[41+w4*ybl[0])/l; 

272 newE3] = (w0*kbl[15]+wl*kbltl3]+w2*kbl[73+w3*kbl[4]'»-w4*kbl[0]) /I; 

273 break; 
274 

275 case 18: 

276 w0-(l-x3); wl-(x3-x2); w2-(x2-x4); w3"(x4-xl) ; w4«xl; 

277 new[03=(wO*cbltl5]+wl*cbl[133+w2*cbl[7]+w3*cbl[13+w4*cbl[0] 

278 new[l]«(w0*mbl[15]+wl*mbl[13]+w2*mbl[7]+w3*rtibl[l]+w4*inbl[0] ) /I; 

279 new[2]-(w0*ybl[15)+wl*ybi[13]-i-w2*ybl[7]+w3*ybl[l]+w4*ybl[0])/l; 

280 new[3]^(w0*kbl[153+wl*kbi[13]+w2*kbl(7]+w3*kbl[l]+w4*kbl[0]) /I; 

281 break; 
282 

283 case 19: 

284 w0={l-x4); wl«{x4-xl); w2«(xl-x2); w3-(x2-x3); w4=x3; 

285 new(0] = (w0*cbl[15]+wl*cbl[14]+w2*cbl(8]+w3*cbl [3] +w4*cbl [0] ) /I; 

286 new[l]«(w0*mbl[15]+wl*nibl[14]+w2*mbl(8]+w3*nibl[3]+w4*m^ 

287 new(2]«(w0*ybiri5]+wl*ybl[14]+w2*ybl(8]+w3*ybl[3]+w4*ybl[0]) /I; 

288 new[33 = (w0*kbl[153-*-wl*kbl[14]+w2*kbl[8]+w3*kbl [3]+w4*kbl [0] ) /I; 

289 break; 
290 

291 case 20: 

292 wO-Cl-x4); wl-(x4-xl); w2=(xl-x3); w3=.(x3-x2) ; w4=x2; 

293 new [ 0 ] = ( wO* cbl 1 1 5 ] + wl * cbl ( 1 4 ] +w2 *cbl [ 8 ] +w3 *cbl [ 2 ] +*w4 * cbl [ 0 ] ) / 1 ; 

294 new[l]-(w0*ir4)lll5]+wl*inbl[14]+w2*mbl(8]+w3*rribl[2]+w4*wbl 

2 95 new{2] = (wO*ybl[15]+wl*ybl[14 3+w2*ybl[8]+w3*ybl[2]+w4*ybl[0])/l; 

296 new[3)-{w0*kbl[15]+wl*kbl[14)+w2*kbl(8]-fw3*kbl[2]-fw4*kbl[0])/l; 

297 break; 
298 

299 case 21: 

300 w0=(l-x4); wl-(x4-x3); w2«(x3-xl); w3=(xl-x2); w4-«x2; 

301 new[0]-(w0*cbl[15]+wl*cbl[14]+w2*cbl[5]+w3*cbl[2]+w4*cbl[0] ) /I; 

302 new[l]-{wO*mbl[15]+wl*nibl[14]+w2*mbl(5]+w3*rTibl[2]+w4*rnbl[0]) /I; 

303 new[2] = (w0*ybl[15]-t-wl*ybltl4]+w2*ybl[5]+w3*yblE2]+w4*ybl(0)) /I; 

304 newt3]-(w0*kbl[15]+wl*kbl[14 3+w2*kbl(5]+w3*kbl[2]+w4*kbl(0))/l; 

305 break; 
306 

307 case 22: 

308 w0=(l-x4); wl«(x4-x3); w2=(x3-x2); w3=>(x2-xl); w4=xl; 

309 new[03-(wO*cbl[151+wl*cbl[14]+w2*cbl[5]+w3*cbl[l]+w4*cbl[0]) /I; 

310 new[13=(w0*inbl[151+wl*mbi[14]+w2*rnbl(5]+w3*mbl[l]+w4*rnbl[0]) /I; 

311 new[2]-(w0*ybl(15]+wl*ybl[14]+w2*ybl[5]+w3*ybl[13+w4*ybl[0)) /I; 

312 newf3]=(w0*kbl(15)+wl*kbl[14]+w2*kbl[5]+w3*kbltl3+w4*kbl[0] ) /I; 

313 break; 
314 

315 case 23: 
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315 wO=(l-x4); wl-(x4-x2); w2-(x2-xl); w3-(xl-x3); w4»x3; 

317 newE0]-(w0*cbl[15]+wl*cbltl4]+w2*cbl(6]+w3*cbll3]+w4*cbl[0])/l; 

318 new[l]-(w0*inbl(a5]+wl*mbltl43+w2*inbl[6]+w3*inbl[3]+w4*ir^^ 

319 new[2] = (w0*ybl[15]+wl*ybl[14]+w2*ybl[63+w3*ybl[3]+w4*ybl[0])/l; 

320 new[3] = (v0*Jcbl[15]+wl*)cbl[14]+w2*Xbl[63+w3*)cbl[3]+w4*3cblt0])/l; 

321 break; 
322 

323 case 24: 

324 wO=(l-x4); wl«=(x4-x2); w2-(x2-x3); w3«(x3-xl); w4=xl; 

325 new[0] = (w0*cbl(15]+wl*cbl[14]+w2*cbl(6]+w3*cbl(l]+w4*cbl[0])/l; 
32 6 new[ll-(w0*inbl[15]+wl*mbltl4]+w2*n4>lt6]+w3*mbl[ll+w4*inbl[0))/l; 

327 new[2] = (wO*ybl(15]+wl*yblC14]+w2*ybl[63+w3*ybl(13+w4*ybl(0])/l; 

328 new[3)«(wO*kblC153+wl*kbl[141+w2*kblt6]+w3*kblll3+w4*kbl[0])/l; 
32 9 break; 

330 

331 } 
332 

333 return; 
334 

335 } 
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What is claimed is: 

1. A method for matching the color input signals of at least 
a first display device and at least a second display device, the 
first and second display devices having an input, compris- 
ing: 

(I) creating a transform LUT for converting input color 
values to output color values by: 

(1) producing a first preselected plurality of input color 
values; 

(2) using said plurality of input values to display a 
plurality of color patches in said first and said second 
display devices; 

(3) obtaining a colorimeiric value of each of the dis- 
played color patches in each of the display devices 
and using said colorimetric values to create a first 
and a second model for the first and second devices 
respectively correlating preselected input color val- 
ues to displayed colorimetric values for each of the 
two devices; 

(4) inputting to said first and second models a second 
preselected plurality of color values; 

(5) comparing an output value signal received from the 
first model to an output value signal received from 
the second model to obtain an error signal indicative 

of the difference between the two output value ^ 
signals from each of said second plurality of color 
values; 

(6) using the error signal to modify the input color 
values to the second model and again comparing the 
output value signal received from the first model to 
the output value signal received from the second 
model to obtain a new error signal; 

(7) repeating the process of steps (5) and (6) above until 
the error signal is a minimum; and 

(8) using the modified color values obtaining at the 35 
minimum error signal to create the transform LUT 
correlating input values to the modified color values; 

(U) displaying a plurality of source color image values on 

said first display device; 
(ni) using the transform LUT to transform the plurality of 

source color image values to corresponding new image 

values before inputting 
said new image values to the input of the second 
display device; and 

(IV) displaying said new color image values on said 
second display device producing a colorimetric match 
with the source color image values displayed on the 
first display device. 

2. The method according to claim 1, wherein in perform- 
ing the step II, first testing every one of the plurality of 
source color image values; and second 

(a) if the source color image values are included in the 
LUT obtaining from the LUX new image values cor- 
responding to the source color image values; and 

(b) if the source color image values arc not included in the 
LUT, then obtaining by interpolation from adjacent 
source image values in the LUT new image values 
corresponding to the source image values. 

3. The method according to claim 1 or 2, wherein between 
steps (I) (7) and (I) (8) there is the additional step of opti- 
mizing the minimum error. 

4. The method according to claim 1 or 2, wherein prior to 
step (I)(5), the gamut of the output color values of the first 
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model is adjusted to be commensurate with the color values 
gamut of the second model. 

5. The method according to claim 2, wherein the source 
color image values are RGB values and the display devices 
arc CRT displays. 

6. The method according to claim 2, wherein the source 
color image values are YMCK values and the display 
devices are 4-color printing devices. 

7. The method according to claim 2. wherein one of the 
display devices is a CRT display and another of the display 
devices is a 4-color printer and wherein there is included in 
the color printer a color converter for converting RGB color 
values to YMCK color values. 

8. The method according to claim 2, wherein one of the 
display devices is a CRT display and another of the display 
devices is a 4-color printer and wherein there is included in 
the CRT a color converter for converting YMCK color 
values to RGB color values. 

9. Apparatus for generating a transform LUT for convert- 
ing input color values to output color values, comprising: 

(1) digital color image values input means for receiving 
digital input color image values; 

(2) an adder device having (i) a first signal input con- 
nected to the input means, (ii) a second, correction 
error, signal input, and (iii) an output; 

(3) a switch circuit having (i) a first signal input connected 
to the adder device output, (ii) a second control signal 
input, (iii) a first output, and (iv) a second output; 

(4) first and second display model means, each for pro- 
ducing output colorimetric tristimulus image values for 
input color values, each of said model means having an 
input and an output, the second model input connected 
to the first output of the switch circuit; 

(5) means, cormected to the output of said first and second 
model means, for comparing the tristimulus values 
output of the first and second display model means and 
for producing an error signal; 

(6) means for receiving and testing the error signal to 
determine if said signal is a minimum error signal, and 
for outputting a correction error signal to the second, 
correction error, signal input of the adder device and a 
control signal to the second control signal input of the 
switch circuit; and 

(7) means also connected to the input means and the 
second output of the switch circuit for generating a 
transform LUT correlating digital input color values to 
modified color values appearing at the second output of 
the switch circuit when the means for testing deter- 
mines that the error signal is a minimum. 

10. The apparatus according to claim 9, further compris- 
ing means for storing the generated LUT. 

11. Tlie apparatus according to claim 9, further compris- 
ing a gamut mapping means for mapping output values from 
the second model into a color gamut commensurate to thai 
of the first model, the gamut mapping means positioned 
between the output of the second model and the means to 
compare. 

12. The apparatus according to claim 9, further compris- 
ing an optimization means in the error testing means for 
determining if the minimum error is an optimum minimum 
error. 



02/24/2004, EAST Version: 1.4.1 



UNITED STATES PATENT AND TRADEMARK OFFICE 

CERTIFICATE OF CORRECTION 



PATENT NO. : 5.579,031 

DATED : NOVEMBER 26, 1996 

INVENTOR(S) : ^^^^ ZHONGJIE LIANG 

It is certified that error appears In the above-indentified patent and that said Letters Patent Is hereby 
corrected as shown below: 

In the ABSTRACT, line 13: change "A output" to -Output-. 
TN THE CLAIMS : 

Claim 1, line 1: change "input signals" to -image values-. 
Claim 1, line 34: change "obtaining" to -obtained-. 



Signed and Sealed this 
Sixteenth Day of S^tember, 1997 



Attest: 




Attesting Officer 



Commissionrr of Patriilx and TriiJemarkx 



02/24/2004, EAST Version: 1.4.1 



